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I.  SUMMARY  OF  WORK  AND  DESCRIPTION  OF  WORK  ACCOMPLISHED 


This  project  has  involved  a  closely  coordinated  theoretical  and  experimental 
investigation  of  the  non-bonded  interaction  between  boron,  aluminum,  and  other  atoms  with 
the  hydrogen  molecule  and  other  partners.  This  work  has  been  carried  out  in  support  of 
work  in  Air  Force  laboratories  on  the  preparation  of  atom-doped  cryogenic  solids,  which  are 
potential  novel  propellant  systems.  Knowledge  of  the  interaction  potentials  is  necessary  to 
predict  the  properties  of  atom-doped  cryogenic  media,  in  particular  solid  hydrogen. 
Moreover,  these  interaction  potentials  must  be  known  for  both  the  ground  and  excited 
electronic  states,  since  the  primary  diagnostic  tool  to  study  doped  solids  in  the  laboratory  is 
electronic  spectroscopy. 

With  the  support  of  this  grant,  we  have  succeeded  in  observing  the  electronic 
fluorescence  spectrum  of  AI-H2  van  der  Waals  complex  involving  excitation  to  several 
different  excited  electronic  states.  This  extends  previous  similar  experimental  work  on  the 
B-H2  complex.  Of  special  interest  is  the  fact  that  we  are  the  first,  and  thus  far  the  only, 
group  to  have  produced  the  B-H2  and  AI-H2  complexes  and  detected  them  by  any  means. 
We  have  also  observed  the  fluorescence  excitation  spectra  of  a  number  of  other  related  van 
der  Waals  complexes,  involving  either  atoms  or  other  binding  partners.  These  data  have 
been  used  to  infer  binding  energetics  and  geometries. 

These  experimental  observations  have  allowed  us  to  elucidate  the  nature  of  the 
ground  and  excited  potential  energy  surfaces,  in  conjunction  with  the  collaborative 
theoretical  effort.  In  particular,  the  global  potential  energy  surfaces  (PESs)  for  the 
interaction  of  a  ground-state  boron  or  aluminum  atom  with  H2  have  been  determined  by  ab 
initio  electronic  structure  calculations.  These  interaction  potentials  have  been  employed  to 
calculate  the  energetics  and  structure  of  the  binary  complexes,  as  well  as  to  investigate  the 
properties  of  A1  and  B  atoms  doped  in  solid  hydrogen.  The  AI-H2  PESs,  as  well  as 
calculated  Al-Al  interactions,  have  been  employed  in  a  new  theory,  describing  interactions 
between  atoms  mediated  by  the  presence  of  spherical  ligands,  to  model  the 
recombination  of  two  A1  atoms  in  the  doped  solid.  An  extension  of  path-integral 
methodology  has  been  completed  to  allow  the  accurate  description  of  extended  systems  in 
which  more  than  one  potential  energy  surface  is  energetically  accessible. 
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The  following  sections  present  in  greater  detail  these  research  accomplishments  and 
relates  our  work  in  the  studies  of  atom-doped  solids.  More  detailed  descriptions  can  be 
found  in  publications  (listed  in  Sec.  Ill)  which  were  prepared  with  the  support  of  this 
grant. 

A.  Electronic  Spectroscopy  of  the  AI-H2  Complex 

We  have  thoroughly  investigated  the  electronic  spectrum  of  the  AI-H2  and  AI-D2 
complexes  and  have  spectroscopically  characterized  the  excited  states  emanating  from  a 
number  of  the  A1  +  H2  asymptotes.^  ^  In  transitions  to  bound  excited  states,  no  resonance 
fluorescence  is  observed,  as  the  excited  states  decay  by  predissociation  to  lower  potential 
energy  surfaces  correlating  with  lower  A1  +  H2  asymptotes.  In  this  case,  the  electronic 
transitions  could  be  observed  by  action  spectroscopy,  in  which  an  emitting,  excited  A1  state, 
which  was  a  product  of  the  predissociation,  was  detected.  An  extensive  series  of  vibronic 
bands  of  both  the  AI-H2  and  AI-D2  complexes  was  observed  in  the  energy  range  of  the  A1 
M  and  4p  atomic  states.  These  bands  displayed  widely  varying  Lorenztian  widths,  in  some 
cases  allowing  the  rotational  stracture  of  the  bands  to  be  displayed.  We  were  able  to 
organize  many  of  the  bands  into  excited-state  van  der  Waals  progressions,  which  were 
interpreted  with  the  help  of  ab  initio  calculations.  Predissociation  occurs  by  coupling  of 
the  bound  Al(3d,  4p)H2  PESs  with  the  repulsive  I'^A^  Al(4.s)H2  PES. 

The  ab  initio  calculations  predicted  the  existence  of  a  strongly  bound,  partial 
charge-transfer  state,  which  correlates  to  Al(3d)  +  H2.  This  state  was  not  observed  in  our 
initial  study^  of  the  A\(3d,  4p)H2  states  since  we  did  not  scan  far  enough  to  the  red.  In 
further  work,^  we  identified  this  state  by  means  of  a  long  excited-state  progression  in  the 
AI-H2/D2  van  der  Waals  stretch  mode.  We  were  able  to  derive  an  effective  one¬ 
dimensional  potential  energy  curve  as  a  function  of  the  van  der  Waals  stretch  coordinate  in 
the  excited  state. 

In  addition,  we  directly  observed"^  the  electronic  transition  to  the  A1(45)H2  (2^Ai) 
state,  which  manifested  itself  as  a  single  broad,  unstmctured  feature  peaking  ~500  cm"l  to 
the  blue  of  the  of  the  A1 4^  25  *-  3p  atomic  line.  This  is  consistent  with  the  repulsive 

nature  of  the  2'^Ai  Al(4i)H2  PES.  We  also  investigated  the  4^  3p  transition  in  the  AlNe, 
AI-N2,  and  AI-CH4  complexes.^  The  excited  Al(45)Ne  state  was  also  found  to  be 
repulsive.  The  45  states  of  AI-N2  and  AI-CH4  are  bound  and  show  vibrational  stracture  to 


\ 
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the  red  of  the  A1 45  <-  3p  ^Pii2  atomic  line,  and  their  observation  required  the  use  of 

fluorescence  depletion  spectroscopy. 

Our  observations  on  the  binary  AI-H2  complex  have  implications  for  the 
interpretation  of  the  spectrum  of  A1  atom  doped  solid  hydrogen.  Krumrine  et  al? 
investigated  the  properties  of  an  atomic  B  impurity  in  solid  P-H2  through  path-integral 
molecular  dynamics  calculations.  They  predicted  the  B  atomic  35  <—  Ip  transition  to  be 
blue  shifted  by  ~6000  cm“l,  or  500  cm-1  for  each  of  the  12  nearest  neighbor  H2  molecules, 
in  reasonable  assignment  with  the  longest-wavelength  feature  in  the  experimental  absorption 
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spectrum  for  B -doped  solid  hydrogen. 
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The  first  absorption  feature  for  Al-doped  solid  hydrogen  appears  at  330  nm. 

considerably  to  the  blue  of  the  longest-wavelength  transition  (45  <-  3p)  in  the  free  atom  at 

394.5  nm.  From  the  observed^  blue  shift  of  the  45  3p  transition  in  the  binary  AI-0H2 

complex  from  that  of  the  free  atom,  we  estimate  that  this  transition  should  appear  in  the 

doped  solid  at  333  nm,  very  close  to  the  observed  wavelength.  The  next  transition  in  atomic 
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A1  is  to  the  3d  state.  Our  experimental  and  theoretical  work  indicates  that  the  PESs 
emanating  from  this  asymptote  are  attractive.  This  suggests  that  the  3d  <-  3p  transition  in 
the  doped  solid  should  be  red  shifted  from  the  wavelength  (308.2  nm)  of  the  free  atom. 
Hence,  the  observed  330  nm  feature  in  the  spectmm  of  Al-doped  solid  hydrogen  can  be 
assigned  to  the  combined  3d,  As  *-  3p  transition.  This  allows  us  to  estimate  the  oscillator 
strength  of  this  transition  as  the  sum  of  the  gas-phase  oscillator  strengths  for  these  two 
transitions,  thus  enabling  the  estimation  of  A1  atomic  concentrations  from  measured 
absorbances.^ 

B.  Potential  Energy  Surfaces  Describing  the  Al(3p)-H2  Interaction 

Because  of  the  orbital  degeneracy  of  the  Al(3/?)  atom,  three  PESs  emanate  from  the 
Al(3/?)  +  H2  asymptote.  Williams  and  Alexander*  °  have  computed  ab  initio  diabatic  PESs 
which  describe  the  interaction  of  the  Al(3p)  atom  with  H2  at  a  fixed  H-H  distance  (r^).  In 
this  effort,  the  computed  adiabatic  PESs  were  transformed  to  an  approximate  diabatic  basis 
for  a  more  convenient  description  of  the  nuclear  dynamics.  Binding  energies  of  the  bend- 
stretch  levels  of  AI-H2/D2  were  also  computed. 

We  have  investigated  in  detail  the  resolved  rotational  structure  displayed  in  some  of 
the  3d  <-  3p  and  Ap  <-  3p  bands  of  AI-H2  and  AI-D2.  Using  ground-state  combination 
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differences  and  employing  optical  optical  double  resonance  fluorescence  depletion 
spectroscopy,  we  determined  the  rotational  constants  and  parity-doubling  parameters  for  the 
ground  bend-stretch  levels  of  Al(3p)-oH2  and  A1(3/j)-/7D2.^  ^  The  derived  spectroscopic 
parameters  were  compared  with  constants  derived  from  fits  to  computed  energies  based  on 
the  ab  initio  PESs  of  Williams  and  Alexander*  **  and  a  new  set  of  PESs  computed  with 
inclusion  of  the  dependence  upon  r,  the  H2  vibrational  coordinate.  The  agreement  with 
experiment  was  excellent  for  the  latter  set  of  PESs,  and  notably  better  than  for  the  fixed-r 
PESs.  Thus,  it  is  necessary  to  take  into  account  the  dependence  of  the  PESs  on  the  diatom 
bond  distance  in  weakly  bound  complexes  of  atoms  with  molecular  hydrogen.  With  the  T- 
shaped  equilibrium  geometry  of  the  ground  electronic  state  of  the  AI-H2  complex,  taking 
into  account  the  vibrational  motion  of  the  H2  molecule  makes  the  vibrationally  averaged 
atom-molecule  separation  smaller,  and  in  closer  agreement  with  experiment.  This  also 
makes  computed  binding  energy  larger  and  hence  presumably  more  accurate. 

C.  Simulations  of  Al  Atom  Doped  Solid  Hydrogen 

Simulations  of  embedded  Al  impurities  in  solid  hydrogen  have  been  carried  out  in 

collaboration  with  the  group  of  Greg  Voth  at  the  University  of  Utah.*  *  In  previous 

theoretical  work  on  embedded  B  atoms  in  solid  hydrogen,  ^  it  was  found  important  to 

incorporate  the  dependence  of  the  interaction  potential  on  the  orientation  of  the  2p  electron. 

The  effect  of  the  electronic  anisotropy  is  even  more  dramatic  for  Al  atom  dopants.  When  a 

nearest-neighbor  vacancy  is  present,  the  Al  atom  becomes  dramatically  shifted  from  the 

original  lattice  site  in  order  to  accommodate  better  the  3p  electron.  ; 

In  this  work  as  well  as  earlier  simulations  of  B-(H2)„  and  B-Ar„  clusters,  the 

nuclear  motion  was  assumed  to  take  place  on  the  lowest  electronically  adiabatic  PES, 
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obtained  by  diagonalization  of  a  sum  of  pairwise  Hamiltonian  matrices.  This  may  be  too 

drastic  an  approximation  for  weakly  bound  complexes  involving  open-shell  atoms,  where 

excited  electronic  states  may  be  accessible  even  at  low  temperature.  We  have  initiated  an 

13 

extension  of  imaginary-time,  path  integral  methods  to  systems  in  which  several  electronic 

14 

states  are  accessed  coherently.  In  a  one-dimensional  model  system  calculation,  we  have 
demonstrated  that  standard  path-integral  techniques  can  be  easily  extended  to  such  systems 
by  exponentiating  the  full  matrix  of  the  potential. 
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D.  Al-Al  Interactions  in  the  Presence  of  Spherical  Ligands 

From  the  viewpoint  of  embedding  either  B  or  A1  in  solid  pH2,  it  is  crucial  to 
consider  how  recombination  of  these  atoms  will  be  mediated  by  the  solvent  molecules. 

Early  work  in  this  direction  by  Voth  and  co-workers^  ^  neglected  the  electronic  anisotropy 
of  both  the  B-H2  and  B-B  interaction.  We  have  recently  shown*  ^  that  it  is  straightforward 
to  go  beyond  this  approximation  by  extending  the  ideas  inherent  in  the  approach  of  Balling 
and  Wright*  ^  to  develop  an  accurate  description  of  the  interaction  of  two  atoms  in  the 
presence  of  multiple  spherical  ligands. 

The  model  is  based  on  accurate  ab  initio  potential  energy  curves  for  all  the  36  AI2 
(or  B2)  electronic  states  which  correlate  with  dissociation  into  ground-state  (^P)  atoms. 

With  a  valence-bond-like  model,  we  transform  the  36  molecular  orbital  states,  which  have  a 
definite  multiplicity  and  Doo/i  spatial  symmetry,  into  a  set  of  36  Cartesian  states  which 
correspond  to  assigning  the  two  3p  electrons  to  approximate  Cartesian  orbitals  centered  on 
either  atom.  It  is  then  easy  to  determine  the  matrix  elements  corresponding  to  the 
interaction  of  each  ^P  atom  with  any  number  of  surrounding  spherical  ligands.  The  lowest 
eigenvalue  of  the  resulting  matrix  defines  the  adiabatic  potential  governing  the  motion  of  the 
atoms. 

In  a  demonstration  calculation,*  ^  we  showed  how  the  interaction  between  the  two  A1 

atoms  is  significantly  modified  by  the  presence  of  the  intervening  pH2  molecules.  At 

moderate  to  long  range  the  interaction  with  the  pH2  molecules  impedes  the  two  3p  electrons 

from  aligning  themselves  correctly  for  the  lowest  electronic  state  of  AI2.  Consequently, 

there  appears  a  significant  barrier  (at  long  range)  to  recombination  of  the  two  A1  atoms.  We 

13. 

are  now  using  this  extended  pairwise  Hamiltonian  model  in  a  full  path  integral  simulation 
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’  of  the  motion  of  two  A1  atoms  embedded  in  solid  pH2. 

E.  Interactions  of  A1  and  H  atoms 

In  view  of  the  interest  generated  by  work  published  in  the  Russian  literature  on 
performance  enhancement  of  propellants  by  the  addition  of  alane  (AIH3)  or  other  aluminum 
hydrides,*^  ^*  we  have  investigatted  bonding  interactions  between  aluminum  and  hydrogen 
atoms.  This  work  has  thus  far  involved  a  study  of  the  metastable  excited  a^II  state  of  AlH 
which  emanates  from  the  ground-state  A1  +  H  dissociation  asymptote.  Previous 
spectroscopic  studies^^  have  precisely  characterized  the  two  singlet  states,  and  A^II, 


emanating  from  the  as5miptote.  We  have  observed  the  b^'L  -  cfiH  electronic  transition  of 

both  AlH  and  AID  in  free-jet  supersonic  beams.  The  d^Ii  state  is  likely  produced  by 

collisional  recombination  of  A1  and  H  atoms  in  the  free-jet  supersonic  expansion. 

We  carried  out  a  spectroscopic  study  of  the  -  a  (0,0)  band  of  both  AlH  and 

AID.  Spectroscopic  constants  describing  the  states  of  both  isotopomers  were  determined. 

We  also  recorded  decay  lifetimes  for  several  rotational/fme-structure  levels  of  the  excited 

state.  These  are  much  shorter  than  in  the  case  of  BH.^^  To  interpret  the  experimental 

observations,  ab  initio  multi-reference,  configuration-interaction  calculations  have  been 

carried  out  for  the  lower  electronic  states  of  AlH.  The  observed  short  excited-state  decay 

lifetimes  are  the  result  of  predissociation  of  the  b^IT  state  by  spin-orbit  interaction  of  this 

state  with  the  repulsive  132+  state  which  arises  from  the  Al(3/7)  +  H(1j)  atomic  asymptote. 

The  excited  state  manifold,  and  the  state  mixing,  is  considerably  more  complex  than  for  the 

24 

excited  states  of  BH,  investigated  earlier  by  Dagdigian,  Yarkony,  and  their  co-workers. 

F.  Complexes  Involving  B  and  Ne  atoms 

25’-28 

In  previous  work  in  our  research  group,  we  investigated  the  evolution  of  the 

2^23^  2^  <-  2s^2p  at  249.8  nm  and  2s2p^  2s^2p  at  208.9  nm  B  atomic 

transitions  in  the  BAr  and  BAr2  complexes.  In  contrast  to  the  corresponding  BAr  state,  the 

BNe  5^2+  state  is  purely  repulsive,  and  the  BNe  5^2+  -  X^II  transition  displays  a  broad, 

asymmetric  feature  peaking  ~270  cm“*  to  the  blue  of  the  B  atomic  line,  characteristic  of  a 

29 

free  •«-  bound  transition. 

We  have  recently  continued  this  work,  successfully  observing  the  3s  2p 
transition  in  the  BNe2  complex.  The  complex  was  prepared  by  a  supersonic  expansion  of 
B  atoms  entrained  in  Ne  at  high  source  backing  pressures  from  a  cooled  beam  source. 

Laser  fluorescence  excitation  revealed  a  second  broad  feature  further  to  the  blue  of  the  BNe 
peak. 

27 

As  in  our  previous  investigation  of  BAr2,  a  variational-Monte-Carlo  (VMC) 
method  was  employed  to  determine  the  wavefunction  and  energy  E  of  the  ground 
vibronic  level  of  BNe2.  A  product  of  Morse  functions  in  the  3  intemuclear  separations  was 
chosen  as  a  variational  approximation  to  the  BNe2  ground  vibronic  wavefunction  since  this 
function  provides  very  good  descriptions  for  the  ground  vibronic  wavefimctions  of  BNe  and 
Ne2.  The  Monte-Carlo  procedure  generates  sets  of  internal  coordinates  which  provide  a 
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representative  sampling  of  the  BNe2  cluster.  We  then  used  these  points  to  simulate, 
semiclassically^^’  the  optical  excitation  spectrum  of  the  complex.  We  found  excellent 
agreement  with  experiment,  both  in  terms  of  the  position  and  the  shape,  of  the  free  ■«-  bound 
absorption  profile. 

G.  Complexes  Involving  Group  14  Atoms 

In  studies  of  non-bonding  interactions  involving  other  atoms,  we  have  explored  non¬ 
bonding  interactions  of  the  Si  atoms  through  laser  spectroscopic  studies  of  the  SiAr  van  der 
Waals  complex.  Specific  impulse  calculations  suggest  that  the  embedding  of  Si  atoms  in 
solid  hydrogen  may  also  offer  advantages  over  /-H2//-O2  as  a  propellant. 

The  electronic  spectrum  of  the  SiAr  complex  in  the  region  of  the  Si  <-  3s^3p'^ 
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3?  transition  near  220.7  nm  was  investigated.  Previous  workers  have  investigated  SiAr 
electronic  states  arising  from  the  excited  valence  Si{3s^3p4s  ^P)  +  At  asymptote.  At  low 
resolution,  a  single  excited-state  (v',0)  progression  of  bands  terminating  in  a  dissociation 
continuum  was  observed.  The  lower  state  of  the  transition  is  the  ground  electronic 
state,  and  the  upper  state  was  assigned  in  the  rotational  analysis  as  32“.  The  ground-state 
dissociation  energy  was  determined  to  be  Dq"  =  178.8  ±  0.4  cm“l.  Ab  initio  calculations  of 
SiAr  and  A3n  electronic  states  correlating  with  the  ground-state  Si(3s^3p‘^  3p)  +  Ar 
asymptote  were  also  carried  out.  Reasonable  agreement  with  the  experimental  results  was 
obtained.  Evidence  for  mixing  of  the  X32“  and  A3n  states  was  obtained  through  the 
determination  of  the  phenomenological  spin-spin  splitting  between  the  Q  =  O'*’  and  1 
components  of  the  X32“  state. 

We  extended  the  study  of  interactions  involving  group.  14  atoms  to  the  Ge  and  Sn 
atoms,^^’  in  order  to  explore  the  role  of  the  large  spin-orbit  interaction  in  moderating  van 
der  Waals  interactions.  ’  Laser  fluorescence  excitation  spectra  of  the  GeAr  complex,  in 

the  vicinity  of  the  individual  fine-stmcture  lines  of  the  Ge  4sMp5s  •«—  4s^4p^  3p  j„ 
atomic  resonance  transition  in  the  region  275  -  260  nm  were  recorded.  Collisional  fine- 
stmcture  relaxation  in  the  free-jet  beam  was  slow,  and  transitions  out  of  all  bound  Ge(4p2 
3f’)Ar  fine-stmcture  levels  were  observed.  The  dissociation  energies  and  vibrationally 
averaged  intemuclear  separations  of  these  states  were  found  to  vary  dramatically  between 
the  different  fine-stmcture  levels,  demonstrating  again  the  importance  of  the  spin-orbit 
interaction  in  determining  the  properties  of  weakly  bound  complexes. 
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Similar  experiments  were  carried  out  on  the  SnAr  complex.  Laser  fluorescence 
excitation  spectra  of  the  SnAr  complex,  in  the  vicinity  of  the  individual  fine-structure  lines 
of  the  Sn  5s^5p6s  jP  *-  5s^5p^  j«  atomic  resonance  transition  in  the  region  317  - 
270  nm  were  recorded.  Transitions  out  of  all  bound  Sn(5p2  3p)Ar  fine-structure  levels 
were  observed.  It  was  not  possible  to  analyze  in  detail  the  rotational  structure  of  the  bands 
since  transitions  from  the  different  SnAr  isotopomers  were  overlapped.  Nevertheless, 
binding  energies  could  be  estimated.  For  both  GeAr  and  SnAr,  the  electronic  states 
associated  with  the  ground  np^  asymptote  are  in  the  case  (c)  limit,  because  of  the  larger 
spin-orbit  interaction.  It  is  found  that  the  large  spin-orbit  interaction  strongly  moderates  the 
strength  of  the  non-bonding  interactions. 

H.  Electronic  Spectrum  of  the  Ga2  Dimer 

We  attempted  also  to  observe  the  electronic  spectmm  of  the  heavier  group  13  Ga 
atom  bound  to  a  hydrogen  molecule.  These  experiments  were  unsuccessful,  but  we  did 
observe  electronic  transitions  in  the  gallium  dimer,  ^  ^  a  species  for  which  no  high- 
resolution  spectroscopic  data  exists.  All  observed  excited  states  are  predissociative  and 
were  detected  by  observation  of  Ga  5s  ->  4p  atomic  emission.  Tentative  electronic  state 
assignments,  based  on  several  ab  initio  studies  of  Ga2,  ’  were  made.  More  calculations 
on  this  molecule  are  needed.  An  upper  limit  for  the  Ga2  dissociation  energy,  D/  <  8891 
cm“l,  was  derived,  which  is  consistent  with  previous  experimental  and  computational 
estimates. 
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